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ABSTRACT. Laplace operators on metric graphs give rise to Klein-Gordon and wave operators. 
Solutions of the Klein-Gordon equation and the wave equation are studied and finite propagation 
speed is established. Massive, free quantum fields are then constructed, whose commutator 
function is just the Klein-Gordon kernel. As a consequence of finite propagation speed Einstein 
causality (local commutativity) holds. Comparison is made with an alternative construction of 
■ free fields involving RT-algebras. 
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1. Introduction 



In the last years the study of quantum systems on networks has received an increasing at- 
tention. They are of interest for possible applications in condensed matter physics. In addition 
interesting mathematical structures appear giving rise to a host of attractive problems, see e.g. 
the articles in |[T3l and further references given there. In this article we study Klein-Gordon 
and wave equations on any metric graph and for any given Laplace operator thereon. We es- 
tablish existence, uniqueness and finite propagation speed for given initial data. In addition we 
J> ; construct free quantum fields on arbitrary metric graphs. The construction of such fields was 

initiated in 01 SI El- The results obtained there were applied to a study of spin transport and 
conductance HEIST]], incorporating additional techniques developed in The main tool 

for the construction of these fields was the use of (a simple version of) RT-algebras IT7l l35ll36ll . 
^ . Also the construction there was limited to relatively simple graphs. The construction we present 

O \ here does not involve RT-algebras and uses only standard and familiar methods of second quan- 

tization. However, we will be able to relate our construction to the RT-construction. Spin will 
not be considered. In order to avoid dealing with infrared problems, we will construct only 
massive and not massless quantum fields. 

We briefly outline our strategy and our results. As a starting point we choose the Hilbert space 
of square integrable functions on the graph as the 1 -particle space. Next we make a choice of 
a self-adjoint Laplacian -A on the graph, which is not necessarily positive. However, -A will 
always be bounded below. To define -A, we follow the discussion in ll26l by specifying bound- 
ary conditions at the vertices of the graph for the operator given as the second derivative acting 
on functions on the graph. Given the Laplacian and a mass m > and motivated by relativistic 
quantum theory, we introduce the energy operator \/-A + m 2 , the d'Alembert operator (wave 
operator) □ = d 2 - A and the Klein-Gordon operator □ + m 2 . Q Unique solutions of the classical 
Klein-Gordon equation for given Cauchy data are then obtained by using 



sin V-A + m 2 t 
V— A + m z 

which is the Klein-Gordon kernel for m > and the wave kernel for m = and which will 
be studied in detail. In particular finite propagation speed will be established. This notion 
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makes sense, since on a metric graph the distance between two points is well defined, so the 
concepts of two events, that is points in space-time, being space-like separated makes sense. 
Finite propagation speed for solutions of the wave equation on smooth manifolds is well studied 
and understood, see e.g (8] [I2j |43] @3). So far for spaces with singularities finite propagation 
speed has been proved only for the case when the singularities are conical j9ll . 

Applying second quantization and a given choice of the Laplacian, we arrive at free fields 
which satisfy the Klein-Gordon equation. They are hermitian as soon as the boundary condi- 
tions defining the Laplacian are chosen to be real, a notion that will be explained below and 
which is equivalent to time reversal invariance in quantum mechanics, when the Laplace opera- 
tor is taken to be a Schrodinger operator. As usual, the non-hermitian scalar fields carry charge. 
For their construction we work with two Laplacians, one for the particle and the other one for 
the antiparticle. They are such that the boundary conditions defining them are the complex con- 
jugates of each other, again a notion that will be explained in due time. Since the commutator 
is actually given by the kernel (11.11 ). Einstein causality (local commutativity) is just another for- 
mulation of finite propagation speed. In other words, we show that the commutator vanishes 
for space-like separated events. Our proof is different from the standard proof of finite prop- 
agation speed on smooth manifolds. Our methods, however, do not allow us to prove finite 
propagation speed and hence Einstein causality in full generality. As a matter of fact, we miss 
those space-like separated events, whose space components are both points in the interior of the 
graph. Theorem [33] gives the precise conditions and statements. 

The article is organized as follows. In Section [2] we summarize several properties of Laplace 
operators on metric graphs in a form needed for the next sections. It includes a detailed discus- 
sion of their (improper) eigenf unctions. In fact since these eigenf unctions give us the integral 
kernel of the Klein-Gordon kernel, some of their properties are crucial for establishing finite 
propagation speed. In addition to recalling several results from Il24ll26*ll28ll30ll . we also estab- 
lish new and relevant ones. This includes the following. Viewing the Laplacian as the Hamil- 
tonian of a quantum dynamical system, there is an associated scattering theory. As it turns out, 
the on-shell scattering matrix enters the eigenfunctions l26ll and hence also the integral kernel 
of the Klein-Gordon kernel. The crucial ingredients in proving finite propagation speed are the 
analytic properties of the S matrix. In the single vertex case the information we gain on the S 
matrix is so detailed, that we are able to establish finite propagation speed even in the case that 
the Laplacian has bound states. 

In Section[3]we discuss classical solutions of the Klein-Gordon and the wave equation. There 
we also formulate the finite propagation speed result, the proof of which is given in Appendix El 
In Section 0] we construct space-time dependent relativistic free fields, both hermitian and non- 
hermitian, that satisfy the Klein-Gordon equation and the same boundary conditions as those 
for the given Laplacian. There we also show, that their commutator function equals (minus) the 
Klein-Gordon kernel (11.11) . The proof of the orthonormality of the improper eigenfunctions of 
the Laplacian is given in Appendix lAl 

2. Laplace Operators on Metric Graphs , their spectral properties and 

their eigenfunctions 

In this section and for the convenience of the reader, we recall the construction of self-adjoint 
Laplace operators on metric graphs in terms of boundary conditions. Also we list several of 
their properties, in particular of their eigenfunctions. They will be needed when we establish 
finite propagation speed and when we construct free fields and discuss some of their properties. 
We start with some elementary concepts from graph theory. The material is mainly taken from 

123 . 

2.1. Basic concepts. 

A finite graph is a 4-tuple Q = (V,X, £, d), where V is a finite set of vertices, X a finite set of 
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internal edges and £ a finite set of external edges. Elements in X U £ are called edges, d is 
a map, which assigns to each internal edge ie Ian ordered pair of (possibly equal) vertices 
d(i):={v\,V2} and to each external edge e e £ a single vertex v. The vertices v\=: d~(i) and 
V2=- d + (i) are called the initial and final vertex of the internal edge i, respectively. The vertex 
v = die) is the initial vertex of the external edge e. If d(i) = {v, v}, that is d~(i) = d + (i), then i is 
called a tadpole. A graph is compact if £ = 0, otherwise it is noncompact. Two vertices v and 
v' are called adjacent if there is an internal edge ie I such that v e d(i) and v' e d(i). A vertex 
■u and the (internal or external) edge j e X U £ are incident if f e <9(j). 

We do not require the map 9 to be injective. In particular, any two vertices are allowed to be 
adjacent to more than one internal edge and two different external edges may be incident with 
the same vertex. If d is injective and d~(i) d + (i) for all i e X, the graph Q is called simple. The 
degree deg(v) of the vertex v is defined as 

deg(v) = |{eG £ | die) = v}\ + \{i e X | = u}| + |{« e J | d + (i) = t>}|, 

that is, it is the number of (internal or external) edges incident with the given vertex v and by 
which every tadpole is counted twice. A vertex is called a boundary vertex if it is incident with 
at least one external edge. The set of all boundary vertices will be denoted by dV such that 
\dV\ < \£\ holds. The vertices not in dV, that is in V; nt = V \ dV are called internal vertices. 

The compact graph Qi nt = (V,X,0,d\j) will be called the interior of the graph Q = (V,X, 
£, d). It is obtained from Q by eliminating all external edges e. Correspondingly, if £ 0, the 
noncompact graph Q ext = (dV, 0, £, c%) is called the exterior of Q. We will view both Q mX and 
<? ext as subgraphs of Q with Q mt n G ext = dV. 

Throughout the whole work we will from now on assume that the graph Q is connected, that 
is, for any v, v' e V there is an ordered sequence {v\ = v, V2, ■ ■ ■ , v n -i,v n = v'} such that any two 
successive vertices in this sequence are adjacent. In particular, this implies that any vertex of 
the graph Q has nonzero degree, that is for any vertex there is at least one edge with which it is 
incident. Q mi is connected if Q is. For connected Q, the graph Q ext is connected if and only if dV 
consists of one vertex only. By definition a single vertex graph is a connected graph which has 
no internal edges, only one vertex, and at least one external edge. The star graph S(v) c £ U X 
associated to the vertex v e V consists of the set of the edges adjacent to v and of the vertex v. 

We will endow the graph with the following metric structure. Any internal edge ie I will be 
associated with an interval = [0, ctj] with ai > such that the initial vertex of i corresponds to 
x = and the final one to x = a^. The open interval If = (0, a^) will be called the interior of the 
edge i. We call the number ai the length of the internal edge i. Any external edge e t £ will 
be associated with a semi-line I e = [0, +oo) whose interior is If = (0, +oo). The set of lengths 
{di}iei, which will also be treated as an element of M) x \, will be denoted by a. A compact or 
noncompact graph Q endowed with a metric structure is called a metric graph (G,q). For the 
purpose of a compact notation we set a e = oo for e e £. The metric structure induces a distance 
function d(p, q) > with the familiar three properties 

• d(p, p) = 

• d(p,q) = d(q,p) 

• d(p, q) < d(p, p') + dip', q) 

for all p,p',q e Q. This defines a topology on iQ,a), such that dip,q) is continuous in both 
variables. For any e, e' e £ we call pdist(e, e') = d(<9(e), die')) the passage distance from the 
external edge I e to the external edge I e >. Thus pdist(e, e') = if and only if die) = die') and 
pdist(e, e') > minj e i ai > 0, whenever 9(e) £ die'), dip, q) > pdist(e, e') holds for any pe I e 
and ge I e > . 
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On the graph Q there is a natural Lebesgue measure dp. In particular there is the Hilbert space 
L 2 (Q) of square integrable functions on Q. We write the scalar product as 



(2.1) {^,4>)g = / ^{p)<t>iP)dp 

Jg 

or simply (tp, (f>), if the context is clear. We write x e Ij = [0, aj] for the coordinate of the point 
p& Q Hp lies on the edge j e £ Ul at the point x and we shall say that the pair (j, x) is the local 
coordinate for p. For short and whenever convenient we will also view (j, x) as a point in Q. A 
complex valued function on the graph, or more precisely on Q \ V, may be considered to be a 
family ^ = {tpj}j e £ui of complex valued functions ipj defined on (0, a,j), so by the convention 
just made tp(j, x) = ipj(x). With this notation the scalar product may be written as 

(ip,(j)}= > / t/jjix^j^dx 







Also we define the derivative ip' = d x ^ of as 

d_ 

dx 

We also introduce the following set of boundary values of if) and its derivative as 

(2.2) ±= {^(0)} ieX , v>'= WMhei 

\{lpi(ai)} ieI J \{-lp'i(ai)}iel; 

The ordering of the set £ is arbitrary but fixed as is the ordering in 1. Given an ordering, in 
(12.21) the boundary values on the external edges come first, then the boundary values at the 
initial vertices and finally the boundary values at the final vertices. Note also that ip' is defined 
in terms of the inward normal derivative, which is intrinsic, that is independent of the special 
choice of the orientation on each of the internal edges. 
The Laplace operator is defined as 

d 2 

(-A AB i;) j (x) = -— ^-(x), je I LI £ 
with boundary conditions 

(2.3) Atp + Bx(/ = 0. 

A and B are {\£\ + 2\T\) x (\£ \ + 2\T\) matrices. For later reference we rewrite this condition as 

(2.4) {A, B) (J^J = 0, 

where (A, B) is the {\£\ + 2\T\) x 2{\£\ + 2|X|) matrix obtained by putting the matrices A and B 
next to each other. So (12.41) is the condition 



(2.5) (Jp\eKer(A,B). 

The operator —Aa,b is self-adjoint if and only if the matrix (^4, B) has maximal rank and 
the matrix AB^ is hermitian. Obviously for any invertible C the pair (CA, CB) gives the 
same boundary conditions since Ker(CA CB) = Ker(A, B). Moreover, with these conditions 
Ker(A, B) is a maximal isotropic subspace M.{A, B) w.r.t. the canonical hermitian symplectic 
form on C 2 ^! 4 " 2 ! 1 !) and all hermitian subspaces can be written in this form, see |[26l . Moreover 
M(A, B) = M(A', B') if and only if A' = CA, B' = CB for some invertible C. For a detailed 
discussion concerning the self-adjointness see |[26l l30l . In addition, if the pair (A, B) satisfies 
these two conditions, so does the complex conjugate pair (A, B) giving rise to the Laplacian 
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A^b a - Let n+(AB^) be the number of positive eigenvalues of AB^ , counting multiplicities. 
The identity 

(2.6) n+{ABt) = n+(AB^) 

is clear. In fact, AB* and AB' actually have the same spectrum. 

Proposition 1. The absolute continuous spectrum of each —Aa,b w the interval [0, oo). It has 
multiplicity equal to the number of external edges, \£\. The number of negative eigenvalues, 
counting multiplicities, is at most n + (AB^) ( < \£ \ + 2|X| ). It is equal to n+(AB^) ifl = 0. 

Below we shall see that the external edges provide a natural labeling for the multiplicities of 
the absolutely continuous spectrum. 

Proof. We claim that all Laplacians —Aa,b are finite rank perturbations of each other, that is 
the difference of two resolvents is always a finite rank operator. To see this, consider the Hilbert 
space 

(2.7) H = H{£,l,a) = H £ ®H x , H e = ® e *eH e , Wz = ©fexWi. 

where H e = L 2 ([0, oo), dx) for all e e £ and Hi = L 2 ([0, oj], dx) for all % e X. Then I?{Q) 3? H. 
By Vj with je £ Ul denote the set of all ipj e TLj such that ipj(x) and its derivative ipj(x) are 
absolutely continuous and ipj( x ) is square integrable. Let V® denote the subset of consisting of 
elements ipj which satisfy 

^-(0) = V/(0) = when j <= £ 

^(0) = ^'(0) = ^j(ai) = 4>'(aj) = when j e X. 

Let A be defined as the second derivative operator, A ^ = V 7 "' with domain 

Z> = © ie eurZ>JcW. 

Then the deficiency index of -A is equal to (|<5| + 2|X|, |£| + 2|X|) and every self-adjoint exten- 
sion is of the form —Aa,b for a suitable boundary condition (A, B). Thus the claim follows by 
general results on self-adjoint extensions, see, e.g., Appendix A in HI and the references quoted 
there. The last statement is just Theorem 3.7 in OTTl . □ 

We elaborate on the sufficient criterion n+iAB^) = for the absence of negative eigenvalues. 
For given boundary condition (^4, B) introduce the meromorphic matrix valued function in k 

(2.8) 6(k; A, B) = -(A + ik5) _1 (^ - ikB). 

Observe that S(k; CA, CB) = 6(k; A, B) holds for all invertible C, so this function depends 
only on the maximal isotropic subspace defined by (A B), S(k; A, B) = 6(k; M{A, B)). 

Lemma 2. ( Il26l , Theorem 2.1;\30\,Theorem 3. 72, EH; Theorem 3.7) (5(k; A, B) exists and is 
unitary for all k > 0. Its poles lie on the imaginary axis. There are no poles on the positive 
imaginary axis if and only if AB^ < and then —Aa,b nas n ° negative eigenvalues. 

The condition A^B < has the following local formulation, see Definition 2.3 in OTTl . in 
terms of vertex quantities and which will be used below. By Proposition 4.2 in [30] for given 
boundary conditions (A, B) there is an invertible C such that the two matrices CA and CB have 
a common block decomposition 

(2.9) CA=($A(v) CB = ($B{v) 

veV veV 

where the pair (A(v), B(v)) gives the boundary conditions at the vertex v. Thus we obtain the 
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Lemma 3. The following block decomposition holds for all k 

(2. 10) 6(k; A, B) = 6(k; B(v)). 

/n particular, if the boundary conditions (A,B) are such that AB^ < 0, then A(v)B(v^ < 
holds for all vertices v and therefore no &(v; k) = (5(k; A(v), B(v)) has poles on the positive 
imaginary axis. 

With the notation just introduced there is the following characterization of k- independence. 

Lemma 4. H24II S(k; A, B) is k-independent if and only if AB^ = and hence if and only if 
A(v)B(v^ = holds for all v e V. 

Alternative characterizations of such boundary conditions are given in (30], Remark 3.9 and 
11241 . proposition 2.4. Thus in the single vertex case all k-independent S-matrices are of the form 

(2.11) S = I-2P 

with P being an orthogonal projector and then S' 1 = = S holds. In combination with theorem 
3.7 in lOTl —Aa,b > follows for such boundary conditions, see also Lemma[2] 

The boundary conditions actually fix the graph. More precisely, given finite intervals I{ (ie T) 
and half lines I e (e e £), and functions tp = {i/jj}jeiue on them, generically the boundary 
condition (12.31 ) given by the pair (^4, B) uniquely fixes the graph Q with a maximal set of vertices, 
such that the boundary conditions are local, see (2611301 for details. 

For given I e £ consider the following solution ^(k) of the stationary Schrodinger equation 
at energy k 2 > 0, 

(2.12) -A AB ^(;k) = kV(;k) 

and of the form 

(e-^dji + SiW)^* forje£ 

(2.13) ^(x;k)=| 

[ aiW)^ + /3(k) jZ e" ikx for j e 2. 

So intuitively we are looking at what happens to an incoming plane wave e~ lkx in channel / 
when it moves through the graph. Observe that choosing the Laplacian -Aa.b as Schrodinger 
operator, quantum mechanically this means that we have free motion away from the vertices. 
The vertices in turn act as beam splitters in a way described by the boundary condition (^4, B). 

The number S(k)ji for j is the transmission amplitude from channel / e £ to channel j e £ 
and S(k)u is the reflection amplitude in channel / e £. So their absolute value squares may 
be interpreted as transmission and reflection probabilities, respectively. The elements S(k)ji 
combine to form the scattering matrix 

S{k) = S AtB (k). 

The "interior" amplitudes a(k)ji=aA,B(^)ji an d fl{k)ji=fiA,B{k)ji are also of interest, since they 
describe how an incoming wave moves through a graph before it is scattered into an outgoing 
channel. 

The condition that ; k) satisfies the boundary condition leads to the solution 

(S{k)\ (I 

(2.14) a(k) \=-Z{kT l {A-ikB) 
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(2.15) Z(k) = Z AB (k) 

X(k) = X(k;a) 



Y(k) = Y(k;a) 
The diagonal \Z\ x \2\ matrices e - are given by 

e ±ika^ = e±ikaj ^ for x 

By construction Z(k; A, B, a) is entire in ke C. For Neumann boundary conditions the scattering 
is trivial, S A =o,B=i(k) = I 

The ip l ( ; k) are not in L 2 (Q) , but rather improper eigenf unctions. Their main properties are 
collected in 

Proposition 5. For fixed k > the ip l (;k) are linearly independent. Any function ip on Q sat- 
isfying — A A ^ B ifj — k 2 ip is a linear combination of these ifj l ( > k), provided k 2 is not a discrete 
eigenvalue of—A AyB . 

The proof will be given in a moment. The next proposition will be play an important role in 
our construction of free quantum fields on the graph Q. Set 

(2.16) = H\ B = {k > | det Z AB (k) = 0}. 

Proposition 6. The improper eigenfunctions ; k) satisfy the the following orthogonality re- 
lations 

(2.17) {i> l (;k),ip l '(;W))=2ir5 lJ/ 5(k-k') k,k'el + \E > . 

For any ke M + \ S > they span the space associated to the absolutely continuous spectrum and 
so the multiplicity of the absolute continuous spectrum equals \£\. In particular, if there are no 
discrete eigenvalues, then the ijr( ; k)form a complete set of improper eigenfunctions of—A A ,B 
in I 2 (C f ). 

That there are no discrete eigenvalues means that i) -A a ,b > 0, ii) there are no positive 
eigenvalues and iii) zero is not an eigenvalue. The proof of (12.171 ) will be given in Appendix 
lAl The remainder follows from the previous proposition. Recalling the notational convention 
(I27TT) . (I27T7T ) reads as 

(2.18) / 4> l (p;k)^ l '(p;k')d P = 2TT8 Ll ,5(k-k'). 

JQ 

For the proof we will need a result concerning the existence of positive (= embedded) eigenval- 
ues. 

Theorem 7. ([26], Theorem 3. 1. 113 11 . Lemma 3.1) —A AtB has a positive eigenvalue E = k 2 if 
and only if k e S > . The multiplicity n(k) is finite. The set S > is discrete and has no finite 
accumulation point in Any eigenfunction to a positive eigenvalue is identically zero on any 
external edge. 

For special boundary conditions, one can obtain many positive eigenvalues, just take for 
example Dirichlet or Neumann boundary conditions everywhere. On the other hand, there are 
also nontrivial boundary conditions, that is ones which do not decouple the external edges from 
the internal ones, and which give positive eigenvalues, see Example 3.2 in [26] and Example 
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4.3 in |29]. Also there are examples with standard boundary conditions (cf. e.g. Example 4.5 
in 11301 for the definition), for which there are positive eigenvalues [23]. 

Corollary 8. The quantities S(k), a(k) and f3(k) depend smoothly on k e K+ \ S > . 

Proof. ZA,B(k) is analytic in ke C, so ZamO^T 1 is smooth in ke R+ \ S > and the claim follows 
from the representation (12.141) . □ 

For further reference we denote by ip^ u for ke E > and 1 < v < n(k) an orthonormal basis of 
the eigenspace with eigenvalue E=k 2 >0. By what has just been proved, each such eigenfunction 
is necessarily of the form 

(2-19) V ; ^( ;E ) = i k,i/ ikx , Kv -ikx f • t 

J [uf e mi + Vj e mx for j e 2. 

The orthonormality condition for fixed k is obviously 



/-o on\ / / k.u ik,u'\ j ( k,u ky , k,u ky 

(2.20) (i/j ' ,tp ' )=6 v y = 2_^)yU i ' u{ Oi + v? at 

+ 2ik ( e2ikai " ( e_2ikai -'))}' 

a quadratic form in the u's and u's. Thus we obtain 

Corollary 9. 77je degeneracy n(k) of any discrete eigenvalue E= k 2 > 0, ?/zatf is k e S > , satisfies 
the bound 

(2.21) re(k) < 2|J|. 
particular S > jj empty when Q is a single vertex graph. 

This result compares with Proposition Q] We turn to a Proof of Proposition [5] Linear inde- 
pendence is clear due to the different occurrence of incoming waves in the different if) ( ; k). 
Assume now that rp satisfies —Aa.b'4' = an d the boundary conditions (!2.3l) .The components 
are necessarily of the form tpj(x) = Uje lkx + Vje~ lkx for all je £Ul. Set 4>=ip-J2ke€ v ki> k { ; k) 
such that also satisfies —Aa,b^ = k 2 4> and the boundary conditions. We have to show that 
= 0. Observe that by construction the components are of the form 

, (x)= /% eik:c > 

3 \<\j(' :kr I <>• ik '\ j'el 

such that ^ contains no incoming waves. Therefore the boundary conditions can be written in 
the form 

fs(k)\ 

(2.22) Z{k) u(k) = 

Wo/ 

with 

S(k) = {s fc } fee£ -, «(k) = {Uj}jeT, V(k) = {vj}jel 

viewed as column vectors. By assumption k e S > , so = u./ = Vj = for all e £ , j e Z, and 
indeed vanishes thus concluding the proof of Proposition [5] 



Theorem 10. ([26] Theorem 3.12, H301I Corollary 3.16) The scattering matrix is unitary for all 
k>0, 

(2.23) S r (fe) t = 5(k)~ 1 . 
addition the identity 

(2.24) 5(-k) = S(k)- 1 
between meromorphic matrix valued functions in k is valid. 
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There are analogous relations for a(k), (3(k) in the form 
Lemma 11. The following identities for meromorphic matrix valued functions in k g C hold 

(2.25) a(-k) = f3(k)S(-k) 

/3(-k) = o(k)5(-k). 

Proof. We will simultaneously also give a new proof of (12.24b - Arrange the components ip l j( ; k) 
as a (\£\ + |Z|) x \£\ matrix ; k), such that the components of ; k) form the I th column. In 
view of (12.131 ). the claims (12.241) and (12.251) combined are equivalent to the relation 

(2.26) ^(;-k) = ^(;k)5(-k) 

as an identity of meromorphic matrix valued functions. Here, by the meromorphic properties of 
S{k), a{k) and (3{k), we view each ; k) as meromorphic in k, that is each component i(^j(x; k) 
is meromorphic in k. So if we define 

(2.27) ^(;k) = V(;-k)S(k) 
we have to show that 

(2.28) ^(;k) = V(;k) 

holds as an identity between meromorphic matrix valued functions. Now -Aa,b^ 1 ( k) = 
k 2 tp l ( ; k) holds. Moreover the boundary values tp\ ; k) and ; k)' of ; k), see (12.2I ). are 
also meromorphic. Since the boundary conditions are satisfied for all k > 0, they also hold for 
all k away from the poles by the identity theorem for analytic functions. Therefore they also 
hold for all ifr( ; -k) and hence also for all %jr( ; k). Similarly -Aa,b^ 1 ( ] k) = k 2 i/; l ( ; k) implies 
-Aa,biP 1 ( ! ~k) = k 2 t/; l ( ; -k) and therefore also -Aa,b4> 1 ( ', k) = k 2 ^ z ( ; k). Again by the iden- 
tity theorem for meromorphic functions it suffices to prove (12.281 ) for all k e 1+ \ S > . But by 
Proposition [5] each if?( ; k) is a linear combination of the tl) k ( ; k). By construction 



(2.29) $(s;k) = 



e- ikx 5 jt + 5(k) jZ e ikx for je£ 

(a(-k)5(k)) ,, e" ika; + (P(-k)S(k)) „ e iWx for j e 1. 



But the eigenfunctions ; k) and tfj l ( ; k) satisfy the same defining properties and so by the 
uniqueness of S(k), a(k) and (3(k) we infer (12.28b - □ 

Remark 12. Since S(k) is meromorphic in k, its unitarity for positive k extends to complex k in 
the form of hermitian analyticity 11111,13911 

(2.30) 5(k) t =5(R) _1 . 
Combined with (12.251 ) this gives 

(2.31) 5(k) f = S(-k). 

In particular S{k) is a hermitian matrix when k is purely imaginary. Since each ;— k) sat- 
isfies — Aa,biP (j — k) = k 2 i() l (;—k) and the boundary conditions (A,B), it has to be a linear 
combination of the ip l ( ; k) and so (12.26!) just provides the explicit form. 

We consider the behavior under complex conjugation. Observe that if (A, B) has maximal 
rank and AB^ is hermitian, then the complex conjugate pair (yi, B) is also of maximal rank and 
AB^ is hermitian. So (A, B) also gives rise to a Laplacian. The following lemma is trivial. 

Lemma 13. 11261 Ifip satisfies the boundary condition (A, B) then the complex conjugate wave 
function ij} satisfies the boundary condition (A, B). 
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In particular, ifip is in the domain of Aa,b> then if) is in the domain of A a B and 



(2-32) -A a ,b^ = -A a ^ 

holds. 

This gives the following nice observation, whose proof we omit. Recall relation (12.61 ) in 
connection with Proposition Q] 

Corollary 14. The spectra of the two Laplacians Aa,b and A^ Q agree. Moreover, if tp is 
an (improper) eigenfunction of—AA,B, then tp is an (improper) eigenfunction of A Ab a for the 
same eigenvalue. 

Let T denote transposition of a matrix. 



Lemma 15. ([29] Theorem 2.2) The following identities between meromorphic matrix valued 
functions hold for arbitrary boundary conditions (A, B) 

(2.33) S^ B {W) = S A , B (W) T 

ot A ^)=JZB~^)S A ,B^) T 
l3 AtB {W) = a A , B ^)S A ,B^) T - 

or equivalently 

(2.34) Va,b( ; k ) = 1>aM ; i<)s A B(kf. 

Proof. We give an alternative proof along the lines used in the proof of LemmaQj] Indeed, with 
the notation used there, define for complex k 



(2.35) V(;k) = VA,s(;k) ^(k)- 1 

where we indicate the dependence on the boundary conditions. The aim is to show 

(2.36) ^(;k) = VAs(;k). 

from which (12.331 ) and (12.341 ) follow. Again by the identity theorem for meromorphic functions 
it suffices to prove this relation for k > 0, for which k 2 is not a discrete eigenvalue. For such 
k by unitarity S^B^) -1 = Sa,b(,^) T an d hence for all k e C, again by the identity theorem. 
By Lemma (fT3l) each %jr( ; k) is an improper eigenfunction of -A{A, B,a) with eigenvalue k 2 
and hence must be a linear combination of the ^ B ( k). By construction, the components of 

tp l ( ; k), k > are of the form 



e- ik ^ 7 + Sa.bQQT 1 ^ for jeS 



(2.37) $(x;k) = 



(aAM^SAM^r 1 )^ + WaA^SaM*)- 1 )^ for j e 1 

't-^5 jl + S A ,B^)i^ x forjeS 

(^a~bW)S a ,b^) T ) e ik " + (Pa~M^)SaM^ T ) e" ikx for j e 1. 

But i/j£ B (;k) and ip l (;k) satisfy the same defining properties and so by the uniqueness of 
S A s (kj, a a 5(k) and @ A s( k ) we infer (I2361) . □ 



By Corollary (fT4]) we know that ip l ( ; k)=ip l A B ( ; k) are eigenf unctions of -A^ B with eigen- 
value k 2 . Relation (12.341 ) tells us that they span the eigenspace of -A Ab for that eigenvalue 



FINITE PROPAGATION SPEED AND FREE QUANTUM FIELDS ON NETWORKS 1 1 

since the ipK B ( ; k) do. We shall make use of this observation when we construct massive, free 
charged fields in Section |4"31 

By definition the boundary conditions given by the pair (A, B) are real if an invertible C 
exists such that the pair (A 1 , B') = (CA, CB) consists of real matrices A 1 and B' . An equivalent 
condition is that there exists an invertible C with CA = A and C'B = B, see ll29l . As a direct 
consequence of Lemmas ITTI and [131 we obtain the following two corollaries 

Corollary 16. For arbitrary boundary conditions (A, B), the relations 

(2.38) a AB (k) = a AB (-k), PaM*) = PaM^ 
hold as identities between matrix valued meromorphic functions in k e C. 
Corollary 17. If the boundary conditions (A, B) are real, then the relations 

(2.39) S(k) = S(-k), /3(k) = a(k)S(-k), o(k) = ^(k)5(-k) 

and hence 

(2.40) a(k) = o(-k), /3(k) = /3(-k) 

are valid as identities between matrix valued meromorphic functions in ke C. 

As a consequence of Lemma [T3l we directly obtain 

Corollary 18. For real boundary conditions (A, B) tp is an eigenfunction of—A.A,B whenever 
tft is. Therefore for a given eigenvalue, the associated eigenspace is spanned by real eigenfunc- 
tions. 

So if for real boundary conditions we choose the eigenfunctions ip ,u to be real, then in the 
notation of (12.191) the relations 

(2.41) U Y=W> keS> 
are valid. Similarly, we can rewrite (12.391 ) as 

Corollary 19. If the boundary conditions are real, then the relation 

(2.42) V(;k) = V>(;k)S'(-k) 

is valid. 

Also (12.301 ) and the first relation in (12.391 ) gives 

Lemma 20. (see [26 ] Corollary 3.2, 1129] Theorem 2.2) If the boundary conditions (A, B) are 
real, then S(k) is a symmetric matrix and so for k purely imaginary the matrix S(k) is real due 
to $23T\) . 

Remark 21. For arbitrary boundary conditions (A, B) the equivalent exponentiated form of 
(12321 is 

(2.43) e iA A , B t^ =e -iAA,B,a*^_ 

If the boundary conditions (A, B) are real and hence Aa,b = ^a b a holds, then (12.431 ) is just the 
statement that time reversal invariance holds. In the single vertex case this invariance combined 
with the hermiticity condition on the field (see below) has been used in [ 5] to prove that 5(k) is 
then a symmetric matrix . 

Combined with (12.1 II ) we obtain 

Corollary 22. For a single vertex graph all k-independent S-matrices resulting from real bound- 
ary conditions are of the form (12.111 ) where P is a real, symmetric and idempotent matrix, P 2 =P. 



12 



R. SCHRADER 



2.2. Negative eigenvalues of the Laplace operator and their eigenfunctions. The operator 
—Aa,b ma y have negative eigenvalues. We introduce the sets 

(2.44) S- = T, A B = {k = ik | k > 0, k 2 = -k 2 is an eigenvalue of - A^} 

S < = B = {k = ik I k>0, k 2 =-K 2 is an eigenvalue of - A^b} 

such that trivially S < c S- and let £ = S- U S > , the set of all discrete eigenvalues. We will 
discuss zero as a possible eigenvalue separately in the next subsection 12.31 Since all Laplace 
operators -Aa,b f° r different (A, B) are finite rank perturbations of each other and since the 
ones with Dirichlet and (or) Neumann boundary conditions are non-negative, S < is a finite set 
and the multiplicity of each eigenvalue is finite. If k 2 = -k 2 < is such an eigenvalue with 
multiplicity n(k), there is a finite, orthonormal basis of eigenfunctions ilr- v , 1 < u < n(k). 
Written in local coordinates they are all necessarily of the form 

The orthonormality condition for fixed ke S < is easily calculated to be 

(2.46) 6 v y = (il}-,^- )=-7^2^ s j S j + 2^{ U i U i a i + v i v i a i 

eeS ieX 

+ m (f u Y ( e2ik " - l ) ( e_2lkai - 0) }• 

In analogy to Corollary [9] we obtain 

Corollary 23. 77ze degeneracy n(k) of any discrete eigenvalue E = k 2 (ke E < ) satisfies the 
bound 

(2.47) n(k) < |£| + 2|Z|. 

After a short calculation, the boundary condition can be brought into the form, compare 
(I2~22l . 

' g k=iK,i^ 

(2.48) Z(k = m) I ;u k=ifi ^ I = 0. 



v 



In case the boundary conditions are real, the ,l> may be chosen to be real, that is the coefficients 
s k ' u , uY and v)' u are all real. 

Recall that there is a canonical Lebesgue measure dp on Q. 5{p, q) is the Dirac 5-function on 
Q with the defining property 



8(p,q)f{q)dq = f(p). 

Remark 24. The arguments may also be reversed to show that £ equals the set of zeros of 
det Z(k) in the set {k e C | Rek = 0, Im k > 0} U R+ and that the k 2 with k e T,form exactly the 
discrete spectrum. As a result there is a completeness relation written as 

1 roc 

(2.49) —J2 dk^(p;k)V'(g;k)+ £ = <5(p, g) P.qeG- 

I J ° ke£,l<i/<n(k) 
77je normalization factor l/2ir is due to (12.181 ). 
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2.2. 1. Bound states and poles of the S matrix in the single vertex case. 

In the single vertex case one can actually say much more about bound states. In fact, we will see 
that they are completely encoded in the S -matrix. Thus the negative eigenvalues are the poles of 
the scattering matrix and the corresponding eigenfunctions may be obtained from the residues 
of the poles. To explain this in detail let Q n denote the single vertex graph with n = \S | half-lines 
meeting at the single vertex v. We will label these half-lines from 1 to n. The scattering matrix 
S(k) now simply equals 6(k) = -{A + ikB^iA - ikB). As shown in EH, see relation (3.23) 
there, the S(k) for different k all commute and as a consequence there is a common spectral 
decomposition 1125) . 

(2.50) 5(k> = s K ( k ) = E ^-^p k = p° - p°° + E P^p k - 

KG J KG j KG Jq 

3 = 3(A, B) is a finite set of different real numbers, including possibly the values n = 0, oo. Also 
3q is the subset, where these elements have been omitted. The P's define a decomposition of 
unity of pairwise orthogonal projectors 

(2.51) E pK = I «x«' {P K ) ] = P K , P K P K ' = P K 5 KK ,. 
kg a 

Thus P K is the orthogonal projection onto the eigenspace of S(k) with eigenvalue (k+iK)/(k-i/c) 
(equal to 1 for k = and equal to -1 for k = oo). The multiplicities are < ns(iK) = trP K . 

S(k) is k-independent if and only if J = {0, oo}, that is J = 0- Then P in (12.111 ) is just P°° 
and P° = I - P°°. Moreover S{k) is invertible if and only if k e D U -i Jo- 
Lemma 25. If the boundary conditions (A, B) are real, then the P K are real, symmetric matri- 
ces. 

Proof. Due to the representation 

P K = lim ^-^S(ir), k*0,oo, 

t^k T + K 

each P K with k g 3q is real and symmetric by Lemma [20] So by the same lemma and with 

r>max KG jK 

T + K 



p o_poc = S{iT) _ y l^ pK 

j-i T-K 

KG Jo 

p0 + pOO = J nxn _ ^ pK 
KG 3o 

are real and symmetric and so are both P° and P°°. □ 

Our next aim is to determine the eigenfunctions out of these data. We will conform to 
our previous notation and show S = T> K = {in \ 0<ks 3o}> suc h that the negative eigenvalues of 
—Aa,b are of the form -k 2 . Given k, there are orthonormal unit vectors s 1K ' u (1 < v < n(ift)) 
in C n which span the eigenspace of P K for the eigenvalue 1 (= Ran P K ), 

(2.52) P s =o K i^s 
and hence 

k + i k 

(2.53) S{k)i K - v = ' 1 " 



k-\K~ 

Observe that the entire set of the s lft '^ is automatically orthonormal by (12.511 ) and (12.521) 



(2.54) »r*r (/ =w^- 
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When the boundary conditions and hence also the projectors are real by the previous lemma, 
these eigenvectors may then be chosen to be real. We define the family of functions ip 1K ' u in 
L 2 {Q n ) in terms of its components as 

(2.55) il>) K,v ix) = sf v V2nc~ KX l<j<n,l<v< n(m). 

The orthonormality of this set follows from the orthonormality (12.541 ) of the s 1K ' u . The next 
result shows how the bound states are encoded in the scattering matrix. 

Proposition 26. Let Q be a single vertex graph. For < k < oo appearing in the spectral 
decomposition (12.501 ) the tp 1K ' v as defined by (12.551 ) are normalized eigenfunctions of —Aa,b 
with eigenvalue —k 2 satisfying the boundary conditions. Conversely, if—K 2 with < k < oo is an 
eigenvalue, then S(k) has a pole at k = in. In particular the multiplicity of each such eigenvalue 
is ns(in) and the number rib of bound states (counting multiplicities) equals 

0<«e 3o 

By this lemma there are at most n = \£\ bound states when Q is a single vertex graph. It is 
easy to construct examples where this upper bound actually is also obtained. 

Proof. As for the boundary values of ip 1 *" 11 and its derivative we have 
and we have to show that 

Ai) iK ' v + B f K - v ' = V2k (A - nB)s lK ' v = 0. 



As established in ESI , see also [30] Proposition 3.7, we may instead of (A, B) equivalently use 
the pair (^(k), B(k)), where 

Mk) = -~ (5(k) - W) . B(k) = ^ (5(k) + W) 

and where k > is arbitrary. Actually by the proof given there, k may be chosen arbitrary in the 
domain of analyticity of S(k) and for which S(k) is invertible. By 

(2.56) det 5(k) = (-l)"s< IO °) J] ( ) 

this is the case if k is chosen outside the set Do U _ i^o- In a moment we shall have occasion to 
make use of this observation. A trivial calculation using (12.53! ) gives 

(A{k) - KB(k))s K ' u = U(-l-^\ 5(k) l n 

~ 2 VV iky k-i« + V ik 

As for the converse let tp # be an eigenfunction of -Aa,biP with eigenvalue —Kq with kq > 0, 
-A J 4_sV ; = ~KQip. Then tp is necessarily of the form i/jj(x) = Cj exp -kqx. Let O^ce C n denote 
the column vector with components Cj. Since tp satisfies the boundary conditions, the relation 

(2.57) (A-n B)c = 
holds or equivalently by the above remarks 

(A(k) - K B(k))c = 0, k £ D U -D 
which when written out gives 

(2.58) 5 (k) C =!^° C . 

k - IKq 
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Thus S(k) has a pole at k = ino g D and P K °c = c For a single vertex graph Z(k) as defined 
by (12.151 ) equals A + ikB. Observe that we used only (12.571 ) to establish (I2.58I ). Therefore the 
condition det (A + iko-B) = for ko ^ 0, cf. Remark|24l is equivalent to ko being a pole for S(k). 
Moreover the discrete spectrum with its multiplicities is given in terms of the scattering matrix 
as 



(2.59) £ = {in | < k g 3o}> n(iK) = ns(iK), < k g 3q. 

As a further, related consequence of this proposition the relation 



□ 



(2.60) P fl = ^ V V 

l<v<.n(\K) 

holds for the matrix elements of P K . 
As for the role of P° and P°° we have 

Lemma 27. The relations 

(2.61) ker^ = RanP°, ker £ = Ran P°° 
hold, so in particular AP° = 0, BP 00 = 0. 

Proof. By our previous discussion ker A = ker A(k) = ker(S'(k) - I) and ker B = ker B{k) = 
ker(5(k) + 1) for k £ O U -D. □ 

The known relation ker ^LLker i? = 0, see Lemma 3.4 in OOl . is of course compatible with 
this result. Consider any piecewise constant function ip, that is a function which is constant on 
each edge. Then ip is completely determined by its boundary values tp. Moreover, if ip g ker A, 
then ijj satisfies the boundary condition (I2.3I ). 

2.3. Zero as an eigenvalue. In this subsection we establish necessary and sufficient conditions 
for -Aa,b to have as an eigenvalue. Let tp ^ be such a square integrable eigenfunction, 
—Aa^b^ = 0- Then necessarily ip e (x) = for all e g £ while tpi(x) = 7« + 5iX for ie I and some 
7i, (5j g C, not all vanishing. So with 7 = 5_= {^ii}i&x e C' 1 ', viewed as column vectors 

and with the notation (12.21 ) 

S - u --(j)- ^(])="(z 

with the diagonal matrix T a = diag{aj}j e j and 

\ ' - 

J tJ 

So the boundary condition (12.41 ) takes the form 

(2.62) (A, B) {^\ = {AUa, BV) (^j = 0. 
Using the condition (12.51 ) and viewing 

Ua 

V 

as a 2(\£\ + 2\T\) x 2|X| matrix we arrive at the 
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Proposition 28. The following relation is valid 

(2.63) dimKer A AB = dim ^Ker(A, B) n Ran (^)) = dimKer(^J7 a , BV). 

Observe that dimKer(A B) = \£\ + 2|J| while dimRan ( .f) = 2|J| since Ker [ T ^ ) = 0. 



Generically subspaces of these dimensions have trivial intersection in a space of dimension 
equal to 2(|£ | + 2|Z|), that is they are transversal . This property remains valid even if one space, 
namely Ker(A, B), is required to be maximal isotropic. As a consequence, for generic boundary 
conditions (A, B) we conclude that —Aa,b does not have zero as an eigenvalue. 

Example 29. Consider the interval [0, a] with Robin boundary conditions at both ends 

(2.64) cos TQtp(0) + sin TQip'(Q) = 0, cos t\ ip(a) - sin T\ip'(a) = 0. 

Then 

(AU a ,BV)=( COSTo Sinr °. 

- \cos T\ a cos T\ — sin T\ / 

has non-trivial kernel if and only if a — tan To — tan n = or cos To = cos T\ = ( Neumann 
boundary conditions). 



2.4. Walk representation of the amplitudes 5(k), a(k) and /3(k). In this section we will pro- 
vide an expansion of the amplitudes S(k), a(k) and (3(k) in terms of walks on the graph. We will 
use this result to give the proof of Theorem [33] in Appendix iBl For the scattering matrix S(k) 
such an expansion was already established in [30]. The extension to a(k) and (3{k) is similar and 
goes as follows. For the convenience of the reader we recall those parts of the notion of a walk 
as introduced in ll30l and extended in OTTl and which are relevant for our purpose. A nontrivial 
walk w on the graph Q from j f £ £ U Z to j e £ U 1 is an ordered sequence formed out of edges 
and vertices 

(2.65) {j,vo,ji,vt, ...,j n ,v n ,j'} 

such that 

(i) ji,---,j n e Z 

(ii) the vertices v e V and v n tV satisfy v e d(j), v e d(ji), v n e d(j'), and v n e d(j n ); 
(hi) for any fee {1, . . . , n — 1} the vertex Vk e V satisfies Vk e d(jk) and vt g d(jk+i); 
(iv) Vk = Vk+\ for some k g {0, . . . , n — 1} if and only if jk is a tadpole. 

When j, j 1 g <S this definition is equivalent to that given in QUI . 
The number n is the combinatorial length |w| com b and the number 

n 

k=l 

is the metric length of the walk w. 

A trivial walk on the graph Q from j' e £ U Z to j e £ U Z is a triple {j, f, /} such that 
?; e and v e 9(j')- Otherwise the walk is called nontrivial. In particular, if d(j) = {vo, v±}, 
then {j,vo,j} and {j,v\,j} are trivial walks, whereas {j, vq,j, vi,j} and {j, v\,j, Vq, j} are 
nontrivial walks of combinatorial length 1 and of metric length aj. Both the combinatorial and 
metric length of a trivial walk are zero. 

We will say that the walk (12.651 ) enters the final edge j through the final vertex vo=vo(w) and 
leaves the initial edge j' through the initial vertex v n = v n (w). A trivial walk {j, v, j'} enters j 
and leaves f through the same vertex v. Assume that the edges j,j' e £ U Z are not tadpoles. 
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The following distance relation holds for a point pinQ with local coordinate (j, x) and the final 
and initial vertices of a walk of the form (12.651) 



(2.66) d(p,v (w)): = 
and similarly 

(2.67) d(q,v n (w)): = \ 



\x if p = {j,x), u (w) = <9 (j), 

\ cij - x if p = (j, x), v (w) = d + {j), 

x' if q^(j',x'), v n (w) = dr(j'), 

ay-x' if q = (f,x% v n M = d + (f). 



The score n(w) of a walk w is the set {rij(w)}j 6 x with rij(w) > being the number of times the 
walk w traverses the internal edge ie I such that 

|w| =^ajni(w) 

holds. Let Wjj', j,j' e £ U X be the (infinite if T ^ 0) set of all walks w on Q from f to j. 
Obviously we have the 

Lemma 30. If p = (j,x) and q = (j',x') and if the edges are not tadpoles, then the 
distance between p and q satisfies 

(2.68) d(p, q) < inf {d{p, ^o(w)) + |w| + d{q, v n {w))) 

with equality if j f. 

Observe that with this notation d(p, t>o(w)) < aj and d(q, v n (-w)) < ay. 
Using relation (3.33) in ll30l . relation (12.141 ) may be rewritten as 

/ S(k) \ fl nxn 

(2.69) a(k) =(I-6(k)T(k)r 1 S(k) mx ,, 
\e- ik ^/3(k)/ \0, 

For the sake of clarity we have indicated the type of matrices with n = \£\,m= \T\. Also 
6(k) = 6(k; A, B), see (EU), and 

T(k) = T(k,a) = 

So S(k) alone is obtained as 




\0 e lK ^ 



Omxn 
Omxn/ 

Alternative ways of obtaining S(k) out of the single vertex scattering matrices ©(k) and the 
metric structure of Q are given in lfTTll22ll29ll4T1l . Analogous relations for the amplitudes a(k) 
and /5(k) are 

(\ 

(2.71) a(k) = (0 mxn l mxm mxm ) (I - 6(k)T(k)r l e(k) mx „ 

\0„ 



/3(k) = (0 mx „ mxm e ik ^) (I - 6(k)T(k))" x S(k) 



nxn 
mxn 



0, 
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As a consequence of relation (12.701) the expansion 

(2.72) 5(k) ee /= J2 5(w;kWe ik l w l 

weW ee , 

with 

(2.73) 5(w;k)« / = n5(Vi;k) ijij _ 1 

l=i 

is valid. S^u; k) is the single vertex scattering matrix obtained from the boundary conditions at 
the vertex v. Also this matrix is indexed by those edges having v in their boundary, that is by 
the edges in the star graph S(v). For this we have to assume that there are no tadpoles, that is 
edges whose endpoints are the same vertex. For the details on the expansion (12.721 ). see QUI . 
But then by the same arguments we also obtain similar expansions for the amplitudes a(k) and 
/3(k). Indeed, for i e Z and e e £ let be the set of walks in Wj e such that v (w) = d^i). W~ 
and W+ are disjoint and W ie = W ie U W+ . Then (1277Tb implies 

(2.74) a{W) ie = V S(w;k) ie e ik H 



weW 



0(k) ie = X! S , (w;k)i e e ikfe+ I W D 

W€>V+ 

with otherwise the same notation as in (I2.73I ). 

3. Classical solutions of the Klein-Gordon and the wave equation 

3.1. Existence and uniqueness of solutions. Fix boundary conditions (A, B) and introduce 
the D'Alembert wave operator 

For given mass m > 0, by definition the Klein-Gordon operator is CU,b + ^ 2 > which we will 
discuss first. 

3.1.1. The Klein-Gordon equation. Our first discussion for the construction of solutions is close 
to the familiar one in the relativistic case. Namely, assume m > to be such that —Aa,b + > 0. 
Then actually there is c>0 such that -A^ s+m 2 >c 2 I holds. Indeed, with e^ ^inf spec-A^B < 
the relation £a,b+^ 2 >0 is valid and so the choice c=\I2{e A,B+m 2 ) does the job. We introduce 
the self-adjoint energy operator 



(3.1) h = h ABm 2 = ^-AA,B + m 2 . 

By what has just been said h > cl, so h has a bounded inverse, < hr 1 < c~ l I. For any / e L 2 {Q) 
define 

(3.2) f (± \p,t) = (e Tiht f)(p) 
which satisfy 

(3.3) ±i^-f i±) (p,t) = hf (± \p,t) 

at 

provided / g V(h). Moreover both f^ ± \p,t) satisfy the Klein-Gordon equation 

(3.4) (n A , B + m 2 ) f {± \p,t) = 

provided the stronger initial condition / e V(-Aa,b) is valid. Indeed, since V(—Aa,b) is 
left invariant under exp(=Fit/i) (h and —Aa,b trivially commute), the functions f (± \p, t) are 
in V{-Aa,b) for all times. Due to the choice of the sign in (13.21) . f^ + \p, t) is called a positive 
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energy solution and /"(p, t) a negative energy solution of the Klein-Gordon equation with 
initial condition / (± '(p, t = 0) = f(p). 

For any gz. L 2 {Q), let g^\p, t) be defined similarly to / (±) (p, t). If in addition ge V{-Aa,b), 
an easy calculation shows that 

(3.5) +i{f ± H;t),d t g {±) (;t))g = (f,g}g 

(3-6) (f (± \;t),d t g T (;t))g = 0. 

holds for all t. In the standard context for the Klein-Gordon equation in Minkowski space 
this result is well known, see e.g. [42j, sec. 3b. In particular, the last relation is read as an 
orthogonality relation between positive and negative energy solutions. 

We can use these observations to solve the initial problem for the hyperbolic differential 
equation defined by the operator Da,b + m 2 within the L 2 context. Indeed, for given /, / with 
/ g V(-Aa,b) = T>(h 2 ) and / e V(h) we will provide a solution f(p, t) to the Klein-Gordon 
equation satisfying the initial conditions 

(3.7) f( p , t = 0) = f(p), d t f(p,t = 0) = f(p). 

Following standard notation, we call the pair (/, /) Cauchy data for the Klein-Gordon equation. 
In fact with the choice 

the function 

(3.8) f(p, t) = {e- iht f i+) )(p) + (e lht f { - ] )(p) 

solves the initial condition (13.71) and satisfies the Klein-Gordon equation. We make the conven- 
tion to say that f(p, t) is a solution for all times if for all t f(,t)e V(—Aa,b) holds, /( , i) is 
twice differentiable w.r.t. t in the strong topology in L 2 (Q) and dtf( , t) e V{h) and finally if 
f(p, t) satisfies the Klein-Gordon equation. Similarly we speak of a solution for small times if 
these properties only hold when \t\ < e for some e > 0. Obviously f(p, t) as given by (13.81 ) is a 
solution for all times. 

In standard contexts there is the well known uniqueness of solutions of hyperbolic differential 
equations for given Cauchy data. The standard proof uses energy conservation, see e.g. Ifl2ll43l 
l44l . In the present context we have 

Proposition 31. For given boundary conditions (A, B) let m>0 be such that —Aa,b + m2 > 0. 
Set h = ^J—Aa,b + m2 and let Cauchy data (/, /) be given with f e T>(— Aa,b) an d f e T>(h). 
Then the solution for small times exists, is unique, therefore extendable to all times and of the 
form (EU). 

Proof. For any solution g(p, t) (for small times) we introduce the energy form 

(3-9) < E(g( , t)) = (d t g( , t), d t g{ ,t))g + (hg( , t), hg( , t))g . 

Since the scalar product ( , )g on L 2 (Q) is positive definite and since h > cl > 0, for given t 
E{g{ , t)) = holds if and only if g{,t) = dtg{ , t) = 0. Also E{g{ , t)) is conserved 

(3. 10) j t E(g( , t)) = (d 2 g{ , t), 8 t g{ ,t))g + (d t g( , t), 8 2 g{ , t))g 

+ (hd t g( , t), hg( , t))g + (hg( , t), hd t g( , t))g 
= -((-A a ,b + m 2 )g{ , t), d t g( , t))g - {8 t g{ , t), (-A a ,b + m 2 )g{ , t))g 

+ (d t g( , t), (-A a ,b + m 2 )g( , t))g + ((-A A , B + m 2 )g{ , t), d t g( , t))g 
= 0. 
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We use this as follows. Let f\{p,t) and /2O?, t) be two solutions for small times for the same 
Cauchy data (/, /) and set g = f\ - f2- By assumption and linearity g(p, t) is also a solution for 
small times. Moreover g has vanishing Cauchy data, g(, t = 0) = dtg( , t = 0) = 0, which implies 
, t = 0)) = 0. But this in turn implies E(g( , t)) = for all small t by (13.101 ) and therefore 
g(,t) = d t g( ,t) = for all small t. □ 

Concerning the existence of solutions for given initial data the positivity condition -Aa.b + 
m 2 > may actually be dropped at the price of stronger domain conditions. To see this, we use 
operator calculus in combination with the spectral theorem to rewrite the solution (13.81 ) to the 
Klein-Gordon equation as 

sin Jit 

(3.11) f (;t) = COS ht f+— f — f 

where both cos ht and sin htlh are bounded self-adjoint operators for all real t. The solutions at 
different times s and t are then related by 

(3-12) /(■,*) = 7 d s f(-,s) 

h 

and we observe that this last relation indeed makes sense without the positivity condition -Aa,b+ 
m 2 > 0. More precisely, for any boundary condition (A, B) and mass m > introduce the Klein- 
Gordon kernel 



sin \/-Aa,b + m 2 t 
(3.13) G A ,B,m 2 (t) = - 



-A a ,b + m 2 

which is well defined by operator calculus. In fact, for fixed t and m > the functions 



sin Vz + m 2 1 , 5 

(3.14) z 1— ► — , 2 1 — > cosvz + m^t 

are entire in 2 e C and bounded and real on the real axis. So both G A ,b ,m 2 W an d <^G A,B,m 2 (t) 
are bounded self-adjoint operators for all t and all m > 0. In order to avoid extra superfluous 
discussion for the case m = we also make the convention 

sin Vk 2 + w- 2 i sin k£ 

(3.15) 



m=0 



\/k 2 + m 2 
To sum up, 

(3-16) K;t) = d t G AtB<w? {t) f + G A , B , m i(t) f 

is well defined for all t. It satisfies the Klein-Gordon equation and solves the initial problem if 
both / and / are in V{-Aa,b)- Also (I3.16l )extends to 

(3.17) f(;t) = G A ,B, m 4t-s)d s f(s), 

valid for all t and s. It generalizes (13 - 1 2b - So far, we have not been able to prove uniqueness of 
the solution in this general case, namely when -A^B+m 2 is not necessarily a positive operator. 

3.1.2. The wave equation. We turn to a discussion of the wave operator CU.b- Consider any 
boundary condition (A B). By the discussion in the previous subsection 

(3-18) /(•, t) = d t G A , B , m ^t) f + G AiW=0 (i) / 

is a solution of the wave equation \Z\A,Bf{p, t) = for given Cauchy data /, / e V{-Aa,b)- 
Concerning uniqueness, there is a result analogous to the one for the Klein-Gordon equation, 
see Proposition [3TJ given as 

Proposition 32. Let the boundary conditions (A, B) be such that —Aa,b is non-negative and 
has no zero eigenvalue. Then the solution (13.181) is the unique solution to the wave equation. 
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In terms of the boundary conditions (^4, B) Proposition Q] gives a sufficient condition for the 
absence of negative eigenvalues, that is n+(AB^) = 0, while Proposition (|28T ) provides necessary 
and sufficient conditions for the absence of zero as an eigenvalue. 

Proof. Again we use the energy function (13- 10b . now with the choice h = ^J-Aa,b > 0. By 
assumption hg = implies g = 0. So again for given t E(g( ,t)) = holds if and only if g(,t) = 
dtg{ ,t) = 0. The proof now proceeds as the one for Proposition [3T1 □ 

3.2. Finite propagation speed. In this subsection we will assume the boundary conditions 
(A, B) to be such that B is empty and that zero is not an eigenvalue of —Aa,b- The aim is to 
analyze support properties of the integral kernel of the operator GA,B,m 2 (t)- When m > we set 
cj(k) = Vk 2 + m 2 . The completeness relation ( 12.49b gives 

1 f°° ,i ,i, \ \ 1 1 / i N sina;(k)t 



1 poo 

(3.19) G AB , m 2(t)(p, q) = -Y J / dk ^(p; k)V z (g; k) 

^ l J-oo 



w(k) 

sinw(k) t 



+ E v> k '»^(<z) ... ■ 

keE,l<^<n(k) UW 

By our convention (13.151) . when m = this simplifies to 

(3.20) G AB ,m^(t)(p, q) = — E J ^ * k)^(g; k)-^- 

k€S,l<^<n(k) 

Observe that due to the self-adjointness of -Aa,b and as is obvious from (13.191 ) and (13.201) . the 
relation 



(3.21) G AtB>m ?(t)(p, q) = G A ,B,m<t)(q<P) 

holds for all m > 0. In addition, due to (12.321) the relation 



(3.22) G AtB>m 2(t)(p, q) = G A B, m 2(t)(q,p) 

is valid. As a consequence, for real boundary conditions G ABtTn 2(t)(q,p) is real. 

We define the space of events to be 1x5 and write an event as (t,p). By definition two events 
(t,p) and (s, q) are space like separated if dip, q)> \ t — s\. 

Theorem 33. Assume one of the following two conditions is satisfied. 

• Q is a single vertex graph (2=0), 

• Q is arbitrary and —A. A- b nas n ° discrete eigenvalues. 

Then for any m > the integral kernel GA t B,m 2 (t ~ S )(P> <Z) vanishes whenever (t,p) and (s, q) 
are space like separated and if in addition at least one of the two points p and q is in Q ex t- 

So far we have not been able to remove the restriction that p or q must lie in in Q ex t- As a 
particular case we obtain 

Corollary 34. Gab,™ 2 (f)(P> <?) vanishes for all p e I e c Q ex t and all q e I e i c Q ex t whenever 
t>0 is smaller than the passage distance, t < pdist(e, e'). 

For the free fields to be constructed in the next section this implies local commutativity (with 
the above restriction). We reformulate finite propagation speed in a more familiar form. For any 
closed subset O of Q and any < d define 

O d = {pe Q\ min dip, q) < d}, 
the closed set of points in Q with distance less or equal to d from O. 
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Corollary 35. Under the conditions of the theorem the following holds for the solution of the 
Klein-Gordon equation (or the wave equation) for given Cauchy data (/, /). 

• If f and f both have support in O c Q sxt , then /(•, t) has support in O'*' for all t. 

• If f and f both have support in O, then supp/(-, t) n Gext c O'*' for all t. 

In particular if both / and / have support on the external edge I e , then /(•, t) vanishes on any 
external edge I e > {e! e) as long as |i| < pdist(e, e'). 

4. Free Quantum Fields on Metric Graphs 

In this section we will construct free fields on the graph Q. The reader is supposed to be 
familiar with the basic concepts of second quantization, see, e.g. lfT9ll20ll42ll46l . Also from 
now on we will assume that the boundary conditions (A, B) are chosen in such a way that there 
are no positive (or zero) eigenvalues of-A^s, that is E > =0 and £=£ < , so bound states are still 
allowed. As a trivial consequence of this assumption, the graph has to have at least one external 
edge, 6 0, since otherwise the entire spectrum is discrete and there are positive eigenvalues. 
Finally we will assume that m > is chosen such that -Aa,b + m 2 > 0. 

4.1. Creation and annihilation operators and the RT-algebra. We introduce the creation 
and annihilation operators 

(4.1) a l (k), a'(k)*, k>0 

a Ku , a Ku \ keS, 1 < v < n(k) 
satisfying the commutation relations 

(4.2) [o , (k),a , '(k')*]=27r<y H /<5(k-k'), [a Ku ,a k ' y *]=5 KW 5 u y 

while all other commutators vanish. These operators act in the bosonic Fock space ^{TL\) with 
7i\ = L 2 (G) as the choice of the 1-particle space, that is 

(4.3) = C®Hi® ■ ■ ■ ®H n ® ■ ■ ■ 

(4.4) H n = H 1 ®H 1 ®H 1 ®H 1 , 



such that 7i n is the n-particle space. ® s denotes the symmetric tensor product. a l (k)* has the 
interpretation of a creation of a particle with wave function tp ( ; k), while a Ku * is the creation 
operator of a particle with (bound state) wave function ijj kiU . The normalization in (14.21) is chosen 
in accordance with (12.171) . ( 12. 19b and (12.461) . For reasons which will become clear in a moment, 
we elaborate on this. By the completeness relation ( 12.491 ) any wave function / e L 2 (Q) has a 
Fourier type expansion of the form 

roc _ _ 

(4.5) /(p) = E / */j(W(p;k)+ E 

with expansion coefficients given as 

(4.6) £ ( k) = -!=/■ v'(p; k)/(p)dp, Hk)=/ Wip)f(p)d P 



'2tT JpeQ Jpsg 

such that the Parseval equality holds in the form 

(4.7) </./>£? = E / dk l/Kk)| 2 + E \f U W 

le£ J ° keS,l<i/<n(k) 



7 

We stick to the standard notational convention in QFT and use to denote the adjoint (only) in this case. 
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holds thus establishing an isometry of Hilbert spaces 

L 2 (Q) « L 2 ([0, oo), dk) ® C Ns 

where 

iV s = n(k) <\£\ + 2\l\ 

keS=S< 

is the total number of bound states, counting multiplicities. With this notation the creation 
operator for a particle with an arbitrary wave function / is of the form 



(4.8) 




and correspondingly its adjoint a(f) is the annihilation operator for the wave function /. 

The (self-adjoint) number operator, the second quantization of the identity operator on the 
one-particle space, is 

(4.9) N = — / dk a'(k)V(k) + a * ak ' U - 

n le£ J ° ke£,l<z/<n(k) 

We define h= ^J-Aa,b + m 2 , see the discussion in Section IXT1 to be the one-particle Hamilton 
operator, so its second quantization is the self-adjoint operator 



(4.10) 




and we observe that u(k) is positive for all k e S < = £ by the choice of m > 0. 
The operator 



(4.11) 




can be given the interpretation of the sum of the absolute value of the momenta of all particles 
in a state of the Fock space which does not contain particles with bound state wave functions. 
Stated more abstractly, let P ac be the orthogonal projector onto the subspace of L 2 {Q) corre- 
sponding to the absolutely continuous spectrum of-A^s- Then P is the second quantization 
of the 1 -particle operator y/-AA,BPac- That there is no proper momentum operator in the fa- 
miliar sense has of course to do with the fact that the configuration space is a graph. So the 
notion of translations in space and with the momentum operator as infinitesimal generator does 
not make sense. But what remains is some kind of absolute value of momentum reminiscent of 
the conservation of the absolute value of the momentum of a (classical) particle under elastic 
scattering. Both N and P commute with Hq and are therefore conserved under time evolution. 

With these preparatory remarks we are now in the position to provide an explicit construction 
of RT (reflection-transmission)-algebras fl7j|35]|3H. The main observation is that k in a\k) and 
a (k)* is positive. So we are free to define creation and annihilation operators also for negative 
k. Indeed, we may set 

(4.12) a'(-k)=^S(kW(k) 

I'eS 

a'(-k)* = £S(-k) z V'(k)*> k>0, 

where we recall the general relation 5(-k) = <S(k) -1 = S(ky valid for all real k ^ 0. With this 
definition the relations (14.121 ) remain valid for k < and then a l (k)* is again the adjoint of a l (k). 
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Since for k > the operator a 1 ' (k)* creates a particle with wave function tp l { ; k), by linearity the 
operator a (— k)* as defined by (14.121 ) creates a particle with wave function 

£S(-kW(p; k), 

l'e£ 

which by (12.261 ) equals ; -k). This gives the first part of the next lemma, while the second 
part follows by an easy calculation. 

Lemma 36. For any k>0 the operator a\— k)* as defined by (14.121) . creates a particle with wave 
function ; — k). The extended family of operators 

{a , (k),o'(k)*} j 

satisfies the commutation relations 

(4.13) [a'(k),/(k'H =5 U '5{k-k!) + S{k) w 5{k + W), -oo<k,k'<oo, 1,1' e S, 
again with all other commutators vanishing. 

Remark 37. This realization of a RT-algebra agrees with the one used in H1 I41II . The construc- 
tion (14.121 ) of the a (— k) and a l (— k)* out of the a 1 (k) and a 1 (k)* is reminiscent of the action 
of the Weyl group in the root space of a Lie algebra, by which any root is obtained from the 
set of positive roots 11141 . A different context, where a (scalar) scattering matrix appears in 
commutation relations, is provided in [ 15 ]. 



le£,-oo<k<oo 



4.2. The free hermitian quantum field. For reasons to become clear in a moment, in this 
subsection the boundary conditions (A, B) will be taken to be real. The field operator, again of 
dimension zero, is defined to be 



(4.14) 



$(tp) = e iH '$(p)e- iH * 



dk 



l 



+ 



2vr VMk) 

E 

keE,l<i/<n(k) 



V> z (p;k)e Mk V(k)* + /i. c.) 



ioj(k)t a W,v * 



+ h.c. 



where h. c. denotes hermitian conjugate. By construction, this field is hermitian and Q(t+s,p) = 
e lH *<3?(s,p)e~ lH * holds. Again we use a similar notational convention as the one used for a local 
description of functions on Q. Thus for its restriction to an edge j and with local coordinate 
(J, x) x g [0, a,j] for a point p there the field is given as 
(4.15) 



Eroo dk 



2tt y / 2w(k) 



e ik ^ jZ + 5(k)^e- 



iWx \ p iw{W)t I 



<!>j(t,x)= < 



, 1 [ c k .^ ikx„kj(k)i n k,u* iUf 



a l (k)* + /i. c. 



36« 



Jo 



oo dk 1 



(k) 77 e" ik:E + /3(k) ?7 e ika;N ) e iw(k)t a'(k)* + /». c. 



J. 



Observe that the ip^' u (p) need not be chosen real. However, the reality of the boundary condi- 
tions comes as follows into play. By Corollary [18] the ip l (p; k) and the i)fi- v {p) are also eigen- 
functions of —Aa,b- Since the boundary conditions are real, we can use Lemma l20l and (14.121 ) 
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to simplify the first terms in (14.151 ) using the RT-algebra notation and m(-k) = lo(K) 
(4.16) 

£ ^ —j== j==== ( (e J % + Whie- iWx ) c*** a\kT + h.c)= 
^ Jo V2vr V2cu(k) \\ J J J J 

[°° — 1 f p i(kx-Mk)f) i fk x* . \ ■ c 

§ r ( (^ e_ikx ' + ^ eikx ') ^ akkr + c - ) = 

£ T "TPt^ (fl^*^ a' (k)* + fc. c. ) , j e J. 
^ V-ao V27T V2tJ(k) v / 

Let O denote the vacuum. 

Proposition 38. The hermitian field 3> satisfies the Klein- Gordon equation 

(n AB +m 2 )<S>(p,t) = 0. 
For all times t the boundary conditions 

A<t(t) + B&(t) = 

are valid in the sense of expectation values in states which are linear combinations of states of 
the form 

i 

with fi e V(-A a ,b + m 2 ). 

Proof. For general boundary conditions (A, B) we recall that if ip satisfies the boundary condi- 
tion (12.31) . then $ satisfies the boundary condition (A, B). As a consequence, if the boundary 
conditions (A, B) are real, then both e^^ip 1 ^ • k) and e lu}( ^ >t ip l ( ; k) satisfy the boundary condi- 
tion (12.31) for all / e £ and all ke R + and the claim follows from the construction of $ and the 
choice of the states. We omit details. □ 

We also introduce the conjugate field 

(4.17) U(p,t) = ^(p,t) = ^-Mp,t). 

at 

Using the completeness relation for the eigenfunctions of —A Ai b in the form ( 12.49b we derive 
the 

Theorem 39. For the boundary conditions (A, B) the equal time commutation relation 

(4.18) [^(p,t),U(q,t)] = i5(p,q), p,q^G 
is valid. 

Observe that this relation fixes the normalization of the field. 

4.3. The free complex quantum field. We now construct a complex field which has the 
advantage of being able to carry (electric) charge. Associated is a particle with that charge and 
an antiparticle with the opposite charge. Accordingly the 1-particle space Hi is chosen to be 
L 2 (Q)@L 2 (Q), the first for a particle and the second for the corresponding antiparticle. 
The 1-particle Hamiltonian h on that space is chosen to be 

(4.19) h = \J-A AB + m 2 © + © ^-A A _ B + m 2 . 

The boundary conditions (A, B) themselves may be chosen arbitrary. To simplify the exposition 
we assume that —A a ,b an d hence also —A^ B has no discrete spectrum, cf. Corollary [14] So in 
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particular -Aa,b > 0, -A^ B > 0. Since the boundary conditions (A B) are not necessarily 
real, relation (12.421 ) need not hold. However, ip( ; k) satisfies the boundary conditions (A, B) by 
Lemma [T3l The creation and annihilation operators for the particles are as before, see (I4.2I ). As 
for the antiparticles, for Ze 6 and k>0 introduce operators b l (k) and their adjoints &'(k)* satisfying 
commutation relations of the same form and commuting with all a 1 (W) and a 1 (k')*. They are 
the annihilation and creation operators for the antiparticle with wave function ip l ( ; k)=t/j l A B ( ; k), 
which we recall differs from ip K ^ ( ; k). Correspondingly we set 

b l (-k)=J2 S(-k) l , l b l '(k)=Y,W)iub l '(k) 

b l (-kr = ]r sRv i^'oo* = E s ^ i ,bl 'w k > 

l'e£ l'e£ 

with S{k) = Sa,b(^)- By (12.261 ) the interpretation is that i/(-k)* creates a particle with wave 
function ; -k). H, the second quantization of h as given by (14.191 ). is 

1 roo 1 poo 

H = — J2 / * W (k) a'(k)V(k) + -E / dk W (k) 6'(k)*6 z (k). 

The field ^ and its adjoint is now given as 
(4.20) 



dk_ 
' x dk 1 



poo AW 1 , 

m,V) = e iH ^(p)e- iH * = £ / Z- /=^ ^ l (P; k) (e iw(k)i 6 ; (k)* + e^ (k)i a'(k) 

rrwo V27T \Jluj{k) \ 



ty.p) = e iH **t(p) e - iH * = V / "L 1 ^(p; k) ( e" Mk) ' &'(k) + e lw(k) * a'(k)* 

^Jo v2vrV 2 ^(k) v 

In local coordinates and in terms of the RT-algebra we can write the field \& (and similarly its 

adjoint) as 

(4.21) 

' * 1 (e iMk) *- ka;) &?'(k)* + e - iMk)t+ka V(k)) , j e £ 



*i(t, x) = 



^/2w(k) 

Use has been made of (12.251) . The motivation for this definition of the one particle Hilbert space 
for the antiparticle, the corresponding 1 -particle Hamiltonian and finally the field \& stems from 

Proposition 40. The field ^(p, t) and its adjoint ^{p, t) satisfy the Klein-Gordon equation 

(\3 + m 2 )^(p,t) = 0, (□ + m 2 )* t (p,t) = 

and the boundary conditions 

AMf) + = = A^(t) + fl*t'(t) 

for all times. 

As in Proposition [38] the last relation holds in the sense of an expectation value in suitable 
states. 

Let C denote charge conjugation, the operation which interchanges particles and antiparticles. 
In addition introduce the antilinear and antiunitary time reversal map T , cf. Remark [21] Then 
there is CT invariance, that is 

(4.22) CT^(p, t)(CT) _1 = *(p, -t) 

holds. 
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4.4. The commutator function. In this subsection we calculate the commutator of the fields. 
For the hermitian field we obtain 

(4.23) 

[mp), Hs, q)] = £ J Q 2^2^ k ^fe k)e Mt - s) - ^(p; W(Q\ k)^^) 

+ E —L- (^■ v {p)^ Ku {q)e luj(m - s) - i) Kv {p)i) k ' v {q)e-^ m - s) ) 
Since the boundary conditions are real, the reality properties 

(4.24) ' ' 

^^( K k^(?;k) = ^^(p;k)^( g ;k); ]T ^(p)^"(?)= £ V k '»^(g) 
Ze£ Ze£ l<i/<n(k) l<i/<n(k) 

hold. Indeed, the first relation is easily derived from (12.42b - To prove the second one, observe 
that for given k e E = S < both sides give the unique integral kernel for the orthogonal projector 
in L 2 {Q) onto the eigenspace of -Aa,b with eigenvalue k 2 < 0. In fact, since the ifr ,u form an 
orthonormal basis in that space, so do their complex conjugates. Inserting the relations (14.241 ) 
into (14.231) gives the first part of the next theorem. The proof of the second part is even easier 
and will therefore be omitted. 

Theorem 41. The commutator for the free hermitian field with real boundary conditions (A, B) 
is given as 

(4.25) [$(t,p),Hs,q)]=-iG AB , m 2(t)(p,q). 

Similarly for the complex field and arbitrary boundary conditions (A, B) the commutators are 

(4.26) [tt(t,p),¥(8,5)t] =-iG AAm 2(t)(p,q), [tf(t, p), q)} = 0. 

The last relation of course also implies ^(t,p)\ \I/(s, qfi = 0. In the Minkowski space 
context it is well known that (up to a sign) the Klein-Gordon kernel equals the commutator 
function, see e.g. fi2l sec. 7c and (14.281) below. So in analogy to the Minkowski space context 
and as a consequence of finite propagation speed we have local commutativity in the form 

Corollary 42. For space like separated events (t, p) and (s, t) the commutators (14.251) and (14.261 ) 

vanish provided as least one of the points p and q lies in Q ex t- 

4.5. Examples. We illustrate our discussion in the context of single vertex graphs with two 
simple examples. First we make the following notational convention. If p has local coordinate 
(i, x) and q the local coordinate (J, y) and for given (A, B) and m we set 

(4.27) Gij(t, x; s, y) = G A B m 2(t - s)(p, q). 

Also we will need the following quantities. Let A(t, x; m) be the usual relativistic commutator 
function of mass m in 1 + 1 space-time dimensions 

(4.28) A(t, x ; m )=^-J ^R) ( e ~ Mk)< - e *) = A<+) (*' x ; m > + A<_> ^' x ' m ) 

(4.29) = riy te sinu,(kK 



oo 



2tt cj(k) 



More explicitly 

(4.30) A(t, x; m) = 



t 2 - x 2 < 

-sign t No(mVt 2 - x 2 ) t 2 - x 2 > 0. 
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iVo is the zero'th Neumann function (a Bessel function of the second kind). For large argument 
it satisfies 

(4.31) N (z) ~ \ — sin(z - vr/4) for 1 < z. 

V irz 

For a more detailed discussion of the commutator function in local coordinates and which will be 
needed in the proof of Theoreml33"1in Appendix iBl introduce the distribution in 0<x, -oo<£<oo 

rf k / e -i^(k)t _ p Mk)t^ k + IK 



1 poo . 

(4.32) W , w >=-/ — e*»(, 



2w(k) V / k-i/t 

= x; m, k) + £> (_) (i, x; m, k) 

00 dk ika . sin u(k)t k + m 



2tt u(k) k-i/c 
with m > and k real, the values k = 0, oo being allowed, that is 

(4.33) T>(t, x; m, k = 0) = A(t, x; m), D(t, x; m, k = oo) = — A(t, x; m). 
By construction D{t, x; m, k) is odd in t. For k oo write 

(4.34) £>(i, x; m, k) = A(t, x; m) + d(t, x; m, k) 

= D {+) (t, x\ m, k) + D ( ~\t, x\ m, k) 
&~'(t, x; m, k) = Ar~\t, x; m) + d^~'(t, x; m, k) 
with the bona fide function 

(4.35) V(t, x; m,K) = — f°° ( - ' 



2m J-oo 2w(k) V / k - w 

= X)^ + \t, x; m, k) + Q^~\t, x; m, k) 

. dk ik sinu;(k)i 1 
= -2ik / — e 



27T LO{k) k-lK 



such that x; m, K) = -ft (+) (-i, x; m, K) = d < - + \t, x; m, k). It is easy to show that for given m 
and k D (±) are uniformly bounded functions of x and t and Holder continuous in both x and t of 
Holder index <1. £> (+) , and £> (_) , are positive and negative energy solutions of the usual 
Klein-Gordon equation respectively 

(d? -d 2 x + m 2 )£ (±) (x, t; m, k) = (d? -d 2 x + m 2 )d {± \x, t; m, k) = 0. 

Moreover, the differential equation 

(4.36) (y~Bx +K ) 0<±> ^' t] m ' ^ = " 2kA(±) ( x > t: m ) 

holds. d^\t, x;m, k) decays at least like |t|~ 1/2 for large t and fixed x. This is well known 
from the theory of Haag-Ruelle scattering theory, see e.g. lfl6l |20l . Sufficient conditions for 
stronger decay are also well known but do not apply here. When < x < \t\ the stationary phase 
approximation gives 

(4.37) D^(t, x;m,K) = (l- v 2 )^ J , e^i^»- 

V / + , . — IK 



'l-v 2 



■ ik ^J2imi\t\ 



with v = xlt and <p(x, t) = mVl - v 2 t. As a function of t and for fixed x the t 112 decay as well 
as the oscillations are visible in numerical computations of £> (±) . 
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Example 43. (The half-line with Robin boundary conditions at the origin) View the positive real 
axis R+ as a single vertex graph with one external edge, \S\=\. All possible boundary conditions 
at the origin giving rise to self-adjoint Copulations are the Robin boundary conditions and 
which are real 

(4.38) cos r ^(0) + sin r ^'(0) = 0, < r < vr. 

They interpolate between Dirichlet ( sin r = 0) and Normans ( cos r = 0) boundary conditions. 
Denote the resulting Laplace operator by — A r . The scattering matrix is now just a function 

cos r — ik sin r 



(4.39) 5 T (k)=- . 

cos r + iKsm r 

satisfying S T (-k) = SV(k)" 1 for k e C and being of modulus 1 for k e R \ {0}, as it should. 
There is a pole of S T (k) at k = i cot r. So for cot r < this pole lies in the lower k-half-plane 
(the second physical sheet). Then there is no bound state and — A T > 0. Conversely cot r > 
gives rise to a pole of S(k) in the upper half-plane at k = i cot r and correspondingly there is one 
bound state with (normalized and real) bound state wave function 



(4.40) ipbAx) = V2cotre- cotra; 
and with bound state energy 

(4.41) e T = -cot 2 r<0. 

As a consequence -A T > Sb- Observe that 5 T (k) is real on the imaginary axis, as should be by 
Remark [121 Note also agreement with Lemma |2] and Proposition Q] In fact in the present case 
AB^ = cos r sin r = cot r sin 2 r. 

By our general discussion the improper eigenfunctions in this example are given as 

(4.42) ip r (x;k) = e- ikx + S T (kW Wx , k>0. 

This set is complete if cot r < while for cot r > this set combined with the bound state wave 
function (14.401 ) forms a complete set. For finite cot r^0 and with the condition m>max(0, cot r), 
such that u(i cot r) > for the mass, we obtain (0 < x, y) 

(4.43) G(t, x; s, y) = —A(t — s,x — y; m) — A(t — s,x + y; m) — t)(i — s,x + y;m, cot r) 

+ 9(oot r)2 cot r ^(icc>tr)-(t- g )) e _ cotr(g+i/) 
uj(i cot r) 

B is the Heaviside step function. 

This example also provides a nice illustration to a long standing problem, namely to what 
extent the scattering matrix is determined by the cross section |[l^[T8l[32l[33l[34l[38ll . Define 
the scattering amplitude T T (k) by S T (k) = 1 + 2i T T (k), that is 

(4.44) T r (k) = -i COS7 



cos r + ik sin r 

The knowledge of |T T (k)| 2 for all k > only fixes sin 2 r. An additional information, namely 
whether there is a bound state or not, is needed to fix r itself. A way to overcome this dilemma 
and to solve this inverse problem in the present context of quantum graphs has been proposed 
in 1271 

The next example is the single vertex graph with two external lines which may also be viewed 
as the real line with the origin as a distinguished point. As boundary conditions we take the the 
one describing the 5-potential of strength A at the origin. This is a very popular model for 
describing a pointlike impurity. 
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Example 44. (The single vertex graph with two external edges (n= \£\ = 2) and with a boundary 
condition describing the 5 — potential on the line) 

The graph is obtained by considering two copies of M+ with their origins identified. The real 
boundary conditions are given as 

The choices A < and A > describe an attractive and a repulsive 5-potential on R respectively. 
The resulting on-shell scattering matrix is a symmetric 2x2 matrix 

1 /-iA 2k \ 2k-iA 1/1 1\ 1/1 -V 



(4.45) 5 A (k) 2k + iA ^ 2k _ iA y 2k + iA 2 Vl V 2 V" 1 1 - 

The second expression gives the spectral decomposition (12.501 ) of the scattering matrix for this 
example, that is P° = and 

(4.46, p-=i(; j). p^i^ -\ 

It has the additional symmetry 

(4-47) S A (k)=(? iWk/° ' 



i o/ AV ' Vi o 



describing invariance of the boundary conditions under the interchange of the two edges. Using 



local coordinates we arrange the components iftAx; k) (j, I = 1, 2; x > 0) of the two improper 



eigenf unctions -0 ( ; k) as a 2 x 2 matrix 



/„-ikx _ iA ika: 2k ika 

L 2k+iA c 2k+iA c 



, 2k p ik:r „-ikx _ iA p ika' 

\ 2k+iA c c 2k+iA c 



Like the S-matrix this matrix is symmetric, reflecting the parity invariance of the 5-potential. 
Also ordinary plane waves appeal - when A = 0, as they should. The relation (12.261) is easily 
verified. In the attractive case A < there is a bound state with bound state energy e\ = -A 2 /4. 
The two local components of the bound state wave function are both of the form 

(4.48) ^( X )=J-^ j = l, 2 . 



AB^ = 



Observe that 

'0 N 
v A y 

and recall again Lemma[2]and Proposition Q] concerning the number of bound states. 

With the notational convention (14.271) Gj^ B rn 2 in local coordinates can be written as a 2 x 2 
matrix in the form, see (1B.2I ). 

(4.49) G(t,x;s,y)=-A(t-s,x-y;m) Q - A(t-s,x + y;m) K J 

-d(t-s,x + y; m, -A/2) j 

™ iv x ,sinM-iA/2)(r-s)) m*m 1 (1 1 

+ G(-A)(-A) — — e 2 - 

w(-iA/2) 2 V 1 1 

Acknowledgments. The author wants to thank L. Faddeev, M. Karowski, V. Kostrykin, and A. 
Sedrakyan for stimulating and helpful comments. 
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Appendix A. Proof of Relation (12.171) 

The general idea of proof follows a familiar route, see, e.g. B31 . However, the boundary con- 
ditions defining the Laplacian enter in a simple but crucial way, which warrant a more detailed 
discussion. In addition the regularity of the scattering matrix S(k) for k > away from S > will 
be used. For given R > 0, let Q R be the set obtained from Q by deleting from any external edge e 
all points with distance larger than R from its initial vertex v e . On each edge e we introduce an 
extra vertex at distance R from v e and denoted by v e>R . Obviously Q R is a compact graph and 
a closed subset of Q. In particular Q R has no external edges and hence is compact. The set of 
vertices of Q R is given as 

VQr ; = VU {v e<R } ee £, 

In other words, Q R is obtained from Q by removing the external edges ee £, each isomorphic to 
the half-line [0, oo), and replacing each of them by a closed interval of the form [0, R], where the 
vertex v e = d(e) corresponds to Oe [0, R] and the new vertex v e>R to -Re [0, R]. Correspondingly 
there is a Hilbert space L 2 {Q R ) with scalar product denoted by ( , ) R . By restriction any function 
/ on Q defines a function on Q R also denoted by /. In this way any element in L 2 (Q) defines an 
element in L 2 {Q R ) and 

lim (f,g) R = (f,g) 

R— »oc 

clearly holds for any /,ge L 2 (Q). As for the claim (12.171 ). the functions ip l ( ; k) are elements in 
each L 2 (Q R ) but not of L 2 (Q), as already mentioned. Now we write 

; k), /( ; k')) R = -£2^72 ((A AB ^( ; k), /( ; k')) R - ; k), A AB /( ; *))r) ■ 

and perform a partial integration. Since the functions ijr( ; k) satisfy the boundary conditions, 
what remains are only contributions from tp l ( ; k) and its first derivative at the vertices v e>R . We 
now observe 

^W; k) = tf, l e (x = R; k) = e- iRk 5 le + S(k) el e iRk 

^\v e , R ; k) = ^{{x = R;k) = -ik e ~ im 5 le + ikS(k) el e iRk 

ax ax 

and obtain 

; k), /( ; k')) R = -^^2 E ((ike iKk 5 /e - ik5(kke" iRk ) (z iRW 5 Ve + S(k') eV e iRk ' 

- (e R % e +Whi^ Rk ) (-ik'e- iRW 5 lle +ik>s(k>) el ,e Rk ' 



k-k' 



e iR{k - k) 5u'-Y,SMeiS(k>) el , e - 

\ ee£ / 

+ kTF (^ e_l/?(k+k,) + S(k') u/ e iR ^ 



Since k + k' > the second term on the r.h.s. vanishes for R — > oo in the sense of distributions 
by the Riemann-Lebesgue lemma. As for the first term write 

(A.1) -^7 {e lRik - k \i' ~ E WhiS(k') el ,e- iR ^) = 
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Here the first term converges in the sense of distributions to 2ir5(k - W)5w as R — > oo. As for 
the second term we use the unitarity of 5(k) to write 

Sw ~ E WhiS(k') el , = E Whi (S(k) el , - S(k') eV ) 

By Corollary [8] all matrix elements of S*(k) are differentiable functions of ke R+\ S > . Since all 
matrix elements also are bounded by 1 due to unitarity, we have the estimate 

\S(k) u > - S(W)u>\ < const- |k-k'|, k,k / GM + \S > , 1,1' e £ 

whenever |k - k'| is small. Observe that M+ \ S > is a union of open, pairwise disjoint intervals. 
This gives the estimate 



l^-E e . g g(k)e^(kQ e H < k(k ' GM+ \ S >, l,l'e£, 

|k-k'| 

again whenever | k - k' | is small. Therefore and again by the Riemann-Lebesgue lemma the 
second term in (IA.1I) tends to zero as R — > cxd. 

Appendix B. Proof of Theorem 

B. 1. Proof of Theorem [33] in the single vertex case. In the single vertex case, besides a proof 
of the theorem, in this appendix we will provide a detailed analysis of the Klein-Gordon kernel 
when written in local coordinates, see the convention (14.27b - We obtain 



J-oo 2ir lo(K) 



(B.l) Gij(t,x;s,y) = -A(t- s,x-y;m 

+ V 2izP K c~ K(x+v) sin M i/c W ~ g)) 

0<K€ 3 V ' 

where we used Corollary [T7] and relations (12.551 ). ( 12.591 ), and (12.60b - Recall also the convention 
(13.151 ) for the case m = 0. We can rewrite this as 



(B.2) Gij(t, x; s, y) = -A(t - s,x-y; m)5ij - A(t - s,x + y;m) (Sij - 2P i 

- 2_j ®tf ^ s, x + y;m, K)Pji 



+ V 2nP K c K(x+v) sin(aj(iK)(t ~ s)) 

n ^ jl w(iK) 

0</te 3o 



Note that P°° may be the zero matrix. We shall use the representation (IB. lb to prove the theo- 
rem. 

For the single vertex graph the distance between two points p and q with local coordinates 
(i, x) and (j, y) is 



(B.3) d(p,q) = d((i,x),(j,y)) = 



\x-y\ 1=3 
x + y i^j- 



As a consequence the first term on the r.h.s. of (IB.ll) vanishes for space-like separations, a 
well known property of the relativistic commutator function. As for the integral in (IB.ll ) insert 
the relation (12.50b . We observe that d((i, x), (j, y)) < x + y is always valid, so for space-like 
separations x+y>\t-s\ holds and thus we can deform the integral from -cxd to +oo to the integral 
from -oo + ip to +oo + ip for arbitrary p > 0. Indeed, by the analyticity of the first function in 
(13.14b we can apply Cauchy's theorem. During this deformation we pick up a residue at each of 
the poles k = in with < k < p. Each such term, however, is compensated by the corresponding 
term in the sum in (IB. lb . When we let p — > +oo, we claim that the integral from -cxd + ip to 
+cxd + ip vanishes. To see this,view the function k i— ► u;(k) as analytic in the cut (open) upper 
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k-half -plane with a cut from im to ioo, see Figure Q] In this cut upper k-half-plane, the estimate 
Imw(k) < Imk holds. 



Imk 



C- 








. V 
im 


J 






Rek 



FlG. 1 . The upper k-half-plane with a cut from im to ioo for the function u>(k). 
C± form the lips of the cut. 



Moreover both functions 
(B.4) 



1 



JwQtitts) 



1 



2iu(k) 



-Mk)(t-s) 



2icj(k) 

are also analytic there and their sum equals 

simV(k)(t - s)) 
^(k) 

there. Furthermore this sum has no discontinuity across the cut, as it should since it is entire 
analytic. Indeed, replace k by the variable m < A<oo via k=iA-e on the left lip C_ and k=iA+e 
on the right lip C+ with e > 0. But on the left lip 



limu;(iA — e) = V X 2 — m 2 

40 



while on the right lip 



lim uj(i\ + e) = — V X 2 — m 2 . 

40 



Using Imw(k) < Im k for k in the upper half plane we can therefore estimate 



(B.5) 



sin(u;(k)(i - s)) 



< 



= Imk|t-s| 

MoT 



in the upper half plane and which combined with 



M*+y) 



_ Q -\mW(x+y) 



proves the claim. This concludes the proof of the Theorem [33] when the graph is a single vertex 
graph. Observe that we have actually proved 



(B.6) 



Gij(t, x; s, y) = -A(t -s,x-y; m)5i 



when x + y > \t — s\, 

(IB .21 ) compares with (IB.6I ). valid when x + y > \t — s\. If at least one of the points p and q is 
far away from the vertex, that is x 1 or y 3> 1, then the last term on the r.h.s. of (IB .21 ) 
becomes exponentially small, uniformly for all times t and s. To sum up, as far as commutators 
are concerned and by comparison with (14.301) . the contribution from U in (IB .21 ) compares with 
the two preceding terms there. 
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Remark 45. We observe from the proof that in the single vertex case the bound state contribu- 
tions in the definition of the fields are necessary in order to obtain locality. A somewhat similar 
observation was made in the context of integrable models in quantum field theory 1211 . There 
it was observed that bound state contributions in the form factors of the Sine-Gordon model 
were crucial for determining the wave-function renormalization constant. Moreover, in the ar- 
ticles El 1401 local commutation relations for certain integrable models were established, see in 
particular relation (54) in EJ, for which also contributions from bound states are relevant. 



B.2. Proof of Theorem [33] for an arbitrary graph when S = 0. We turn to the case of an 
arbitrary graph with the spectral assumption £ = for the Laplacian —Aa,b, that is with the 
assumption that there are no bound states. In local coordinates 
(B.7) 



GiAt, x; s,y)= < 



-A(t -s,x-y; m)<% + JZ §S(k) lj& 



,ik(x+-y) sinMk)(t-s)) 



for i, j g £ 



sinMk)(t-s)) 

w(k) 



for 



el,je £ 



Jky sin(o;(k)(t-s)) 
w(k) 

i,j g J. 



for 



Relation (12.231 ) has been used for the case i, j g £, Corollary [TT]for the case i g 1, j g £. Lemma 
[TT]and Corollary [l6]have been used for the case i, j g 1. Consider first the case j, I g £. The first 
term, the relativistic commutator function, has already been dealt with and vanishes for space- 
like separations. As for the integral we insert the path space expansion (12.721 ) for the scattering 
matrix to obtain the representation 



(B.8) 



E 

we W, 



2^ 



1-3 



ik(x+y+|w|)£^H^ 

(W-s)) 

w(k) 



Here and in what follows we will freely interchange summation and integration. This is permit- 
ted as can be shown with help of Proposition 5.6 in [30] and where one lets the lengths a-i of the 
internal edges become complex with a positive imaginary part. We omit details. 

For events, which are space-like separated, x+y+\ w| > d((i, x), (j, y))> \t—s\ is valid for any 
w g Wij whenever i, j g £. Also by Lemma[2j the assumption S = implies AB^ < which 
in turn implies that A(v)B(v)^ < holds for all vertices v by Lemma [3] This in turn implies 
that each S(v; k), which is of the form -(A(v ) + ikB(v))~ 1 (A(v) - ikB(v)), has no poles and and 
hence is analytic in the upper half plane and polynomially bounded there, again by Lemma [2l 
As a consequence each S(w; k)ji is analytic in the upper half-plane and polynomially bounded. 
These considerations again allow us to make a deformation of the integration over k in (IB.8I ) 
from the real axis (-oo, +oo) to the parallel line (-00 + ip, +oo + ip). Combining the estimate 
(1531 ) with 

e ik(z+</+|w|) -g-ImkCfx+jz+lwl) 



and the polynomial bound of each S(w; k)ji in the limit p — > +oo we obtain a vanishing contri- 
bution. In other words, each summand in (IB.8I ) vanishes. This concludes our discussion of the 
case i, j g £. 
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We turn to the case i e X and j e £ and discuss the integral involving the a(k) and /3(k) 
amplitudes separately. By the walk expansion (12.741) 

(B.9) 

* a(k) . c ik(^) Sin(a;(k)(t- S )) = y r ^ (w . k) . c ik(x+|w| +?; ) sin(^(k)(t - g)) 
i-oo 2vr lJ w(k) w ^,_ 7-oo 2vr ' * J w(k) 

and we observe that d({j, x), (I, y)) < x + |w| + y, holds for all w e VV^, see (12.66b - Hence 
for space-like separation and for each summand we can again deform the integration contour 
to (-00 + ip, +oo + ip) and thus this expression then vanishes when p — > oo. As for the term 
containing the amplitude /3(k), again the walk expansion gives 

(B.10) 

f°° —B(k) c ik(-z +?/ ) SmMk)(t-3)) = y f°° * 5(w . k) . c ik(a ?: -x + |w| +w ) SM^(k)(t-s)) 
J-oo 2?r y w(k) w ^, + i-oo 2?r ' y cj(k) 

ij 

Now x), (l,y)) < cii - x + \ w\ + y holds for all w e W^, cf. again (12.661) . and the previous 
arguments can again be applied. 

In the case j, I e 1, the arguments just used do not work. This is the reason why we have 
been unable to establish finite propagation speed inside the graph, that is in Q m . Indeed, now 
the contour deformation into the upper k-half plane can not be carried out, since f3^ g(-k) will 
have poles in the upper half-plane. Also the walk representation of ^j(-k) for k > does not 
have the form needed to invoke the arguments we have used so far. 

Remark 46. The reason we had to impose the condition 12 = for a general graph is that in the 
presence of bound states we do not (yet) have sufficient control over the matrix valued functions 
S(k), a(k) and (3(fc) at the poles. Recall that in the single vertex case we had Proposition \26\at 
our disposal. However, we expect Einstein causality still to be valid without this condition. 
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